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ABSTRACT 


A mathematical model has been developed to describe the field of vision 
of a pilot seated in aa aircraft. Given the position and orientation of the 
aircraft, along with the geometrical configuration of its windows, and the 
location of an object, the model determines whether the object would be 
v>ithin the pilot's external vision envelope provided by the aircraft's 
windows. The computer program using this model has already been 
implemented and is described herein. 


1 . 


INTRODUCTION 


The National Aeronautics and Space Administration, Wallops Flight 
Center, Wallops Island, Virginia, is conducting studies of the existing 
air traffic patterns at uncontrolled airports. Data acquired are being 
Used to evaluate new air-traffic pattern concepts which may reduce the 
potential for mid-air collision hazard near such airports. The Air 
Traffic Group of Computer Sciences Corporation, Wallops Island, 

Virginia, is providing support to the project by developing computer 
programs for data reduction, developing a data base, and providing 
statistical and mathematical modeling of the data. Approximately three 
thousand radar tracks were acquired at three uncontrolled and three 
controlled airports. A data base called Integrated Data Store (IDS) 
has been developed to store position-velocity time histories of the 
tracks, type of aircraft, airport, runway used, wind speed and direction, 
cloud ceiling, visibility, barometric pressure and other operator comments. 

This report describes the math model developed for analyzing the 
pilot's field of vision from the cockpit of his aircraft. It has been 
noted that the pilot has a very limited field of vision available to him 
because he cannot see straight up, straight back, or straight down; and 
also his ability to detect other aircraft in his vicinity depends upon 
the orientation of his own aircraft. Most mid-air collisions in the vicinity 
of airports occur because either the pilot fails to look or cannot see due 
to his visual restrict^ons. This math model has been developed to 
determine whether or not some exterior point is potentially visible to a 
pilot seated in an aircraft. It should be mentioned that L. S. Joel, 



I 


W. A. Steele, J. J. Fllllben and G. B. Hare have used slmllai* Ideas regarding 
the geometry of the vision envelopt; in their math model.* 

The computer program simulating the model is Initialized by calling 
Subroutine INIT, which loads the common data areas with the window outlines 
for the chosen aircraft. Window outlines were obtained from binocular 
photogiaphs supplied by Mr. P. M. Rich of the Federal Aviation Administra- 
tion's National Flight Evaluation Center at Atlantic City, New Jersey. 

The Subroutine VISION is then called to test the visibility of a particular 
target from a viewing aircraft . The program presently contains the window 
outlines for the following types of aircraft. 


AIRCRAFT TYPE 

MOONEY-21 1 

CHEROKEE-140 2 

CESSNA-172 3 

CESSNA-210 4 

PIPER AZTEC 5 

PIPER 140B 6 

BEECH BARON 7 

AERO-COMMANDER 680E 8 

BOEING 707 9 

DOUGLAS DCS 10 


Space has been reserved in the common area for ten more aircraft types . 


*Natlonal Bureau of Standards Report 1052, entitled "Approaches to Evaluating 
the Effects of VFR Towers on Flow and Safety at Airports." 


2 . 


MODEL FORMULATION 


Since It is the pilot who is at the control, the model equations are 
written in the coordinate system fixed with the aircraft (A). In essence, 
the problem is to locate the point of intersection (I) of the line 
joining the pilot’s eyes (P) to the object (0) and the plane of the window. 

If this point lies in front of the pilot's eyes and falls within the boundary 
of the window, the object should be visible through that particular window; 
otherwise not. 

2.1 COORDINATE SYSTEMS 

Two distinct coordinate systems are used. Figure 1 illustrates the 
system fixed to the ground (G-System) with the origin at any suitable point. 
Figure 2 shows the system fixed to the viewing aircraft (A-Syste ) with the 
origin at any suitable point. 



R, 6, and a are range, elevation, and azimuth of the aircraft. The 


cartesian components are: 



X =■ R cos 6 sin a 


y “ R cos 3 cos a 
z = R sin 3, 

where -90“ £ S 90“ , 

-180“ (0“) < a < 180“ (360“) . 


( 1 ) 



Figure 2. AIRCRAFT COORDINATE SYSTEM (A-SYSTEM) 


r , X, and y are spherical-polar coordinates of a point in A-System. 

The cartesian components of which are written as in (1) with (R, 3,ct) 
replaced by (r,X, y) , respectively. 

Since the locations of both the pilot and the object are 'read' by the 
Same device on the ground, a coordinate transformation is required to locate 
the object in the A-Systera. Figure 3 illustrates the trajectories of the 
aircraft and the object. It may be pointed out that the object is treated 
ttS pd-uc moS s and may be another aircraft or a fixed structure, such as 
a build’ ig or a mountain peak. 


Z(UP) 



Figure 3 


Let (xq, Yo, Zq) are the coordinates of A and (xg, y^, Zg) of 0 in 
G-System. Also, assume that (Xg, y^, Zg) are cuordiuaLes of 0 in A-System. 
The three vectors are related by a coordinate transformation at any given 
time: 


where 
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(T, 0, (j)) are Euler angles between G- and A-Syslems defined by 


(3) 


(A) 


three successive rotations; about Z-axis by 4* (heading), then about new 
X-axis by 0 (pitch), and finally about new Y-axis by (|) "(bank). 



For the sake of generality, let (Xp, Vp, Zp) be the coordinates of the 
pilot's eyes (!•) in the A-System. In the A-System the windows have fixed 
positions. Therefore, the coordinates of each comer of a window are known. 

2.2 SOME GEOMETRICAL FORMULAE 

The coordinates of the point of intersection (I), of the line from 
P to 0 and the plane of the window, will be determined. 


As it may be found in any textbook on solid geometry, the equation of 
a (straight) line passing through two points (x^, y-^, z^) and (X2, V2y 


x-xi y-Yi 

* 2-^1 ” ^ 2-^1 

where Q is a constant. 


Z-Zi 


Zo-Z 


= Q. 


( 5 ) 
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The distance (D) between two points (xj^, y2^, zj^) and (x2» y2* Z2) is given by 

D = Y(x2-Xi) 2 + (y2~yi)2 + (z2-Zj)2. ( 5 ) 


Any linear equation in (x, y, z) defines a plane. The equation ». a 
plane, whose three non-collinear points are (xq^, y^, Zj^), (x2, V2t ^2^ 
(X3, y^, Z3), is given by determinant equation: 


X y z I 

xi y^ z^ 1 

X2 y2 Z2 1 

X3 y3 Z3 


( 7 ) 


1 


The intersection point (x^, y^, of the line defined by ( 5 ) and the 
plane defined by ( 7 ) is given by 



= Q(xa-Xp) 

+ 

Xp. 



“ Q(ya-yp) 

+ 

yp. 

(8) 

"1 

= Q(VV 

+ 

"p* 



where d-axp-by_-cz_ 

Q = , 

a(x-Xp) + b(y-yp) + c(z-Zp) 


with a= yi(z2-Z3) “ Zi(y2"y3) + (y2^3"22y3), 

b= -X]^(z 2 ~Z 3 ) + zi(x 2 -X 3 ) - (^- 2 ^ 3 " 22 * 3 ^ 
c= X]^(y2-y3) - yi(x2-X3) + 

d= X3(y2Z3~Z2y3) - yi(x2Z3-z2X3) + zj^(x2y3-y2X3) . 

The area ( A ) of a triangle, whose comers are located at (x^, y^ , z^) , 


(X2, y2» Z2) and (X3, y3* Z3) , is given by 

U 23 = 1/2 -^k-^ + A2^ + , 

where A3 = X3(y2~y3) - yi(x2-X3) + (x2y3-y2X3), 

A2 = X3<Z2-Z3) - Z 3 (X 2 ~X 3 ) + (X2Z3-Z2X3), 

A 3 = yi(z2-z3) " zi(y2-y3) + (y223-22y3)* 


( 9 ) 


2 . 3 GEOMETRY OF VISION 

To find out whether the intersection point (I) lies within the boundary 
of a window the procedure is described as follows: 



Figure 4 shows a typical window for the sake of illustration. 



Numerals (1, 2, 3, 4, 5) index comers and I is the point of intersection. 
The window and the point I lie in one plane. We form all possible combinations 
out of corners taking two adjacent ones at a time, i.e., (1,2), (2,3), (3,4), 
(4,5), and (5,1). We pick up one pair at a time and the point I each time to 
form as many triangles as there are pairs, i.e., (1,2,1), (2,3,1), (3,4,1), 
( 4 , 5 , 1 ) and (5,1,1). We sum the areas (say, A^) the triangles above. 

Then we calculate the true area (say, A*^) of the window, which is the sum 
of the areas of triangles whose vertices are at (1,2,3), (1,3,4), and (1,4,5). 

It may be noted that, by definition, the innter angle between any two 

adjacent sides of a window may not exceed 180°. Next we evaluate the 
distance (say, D ) between che pilot's eyes (P) and the object (0), and 
also the distance (say, between the point of Intersection (I) and the 

object (0). Quantities A-*- and A^ are calculated using expression (9), while 
and are evaluated using (6). 

2.4 VISIBILITY CRITERIA 

Define a small arbitrary number (say, e) which will provide the model 
a certain tolerance in ‘•he accuracy of the input data and the computer. 

Tlien, the table in Figure 5 establishes the criteria that the object (0) 
would be visible to the pilot (P) through the particular window whose plane 
has the point of intersection (I) . 
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1 

Ijto < pPO 

1 

VISIBLE 

NOT VISIBLE 

d“ > 

NOT VISIBLE 

NOT VISIBLE 

1 

! 




Figure 5. VISIBILITY TABLE 






3. 


NUMERICAL EXAMPLE 


For the purpose of illustration, we have chosen a CESSNA-172, the vision 
envelope of which, as seen by the pilot, is depicted in Figure 6. As it i. ay 
be noted, the windows have shapes complicated enough to be described by 
simple geometry. Therefore, the windows are approximated to simple 
geometric shapes for mathematical convenience, as shown in Figure 7. 

We will assume that the origin of the A-System is located at the 
pilot’s eyes; therefore, 

Xp=yp = zp = 0* 

The position coordinates of the aircraft and the object are given in a 
coordinate system fixed to the runway, which we will choose to be our ground- 
based G-System. The orientation of the aircraft (i.e., pitch, heading, and 
bank angles) are given in the radai -lelaLive coordinate system. Since X-'i 
planes of both the ground and radar are the same, only the heading is 
transformed to the G-System. 

We will now show, using our model, whether an object is visible through 
the segment of the front window, marked (X) in Figure 7. The coordinates of 
the vertices of this window are calculated using (1) in the A-System, and are 
listed in Figure 8. r is of unit length, which seems to be a valid 


( 10 ) 


assumption. 




Figure 6. Vision Envelope for Single-Engine, High -Wing 


1966 CESSNA 172 
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Comer No. (in degrees) (In degrees) • 


-2 

0 

0.000 

0.999 

-0.0348 

-2 

46 

0.718 

0.694 

-0.034 

46 

60 

0.6015 

0.347 

0.719 

65 

0 

0.000 

0.422 

0.912 


Figure 8 


The coordinate of the aircraft and the object in G-System are given as 


yo. = (-2142. 6498, 770) ft, (11) 

(Xg, yg, Zg) = (-11050, 4000, 1200) ft (12) 

The orientation of the aircraft in radar-relative system is 

(4', 0,4>) = (-69°, - 2 \ -19°), (13) 

and in the G-System (runway) is 

(’f, 0, <!>) = (-69°-46.18°, -2°, -19°) 

= (-115°. 18, -20°, -19°), (14) 


where we have corrected for the runway azimuth (46°. 18) as measured by radar. 
The coordinates of the object in the A-System are evaluated using (3) and (4) 
With values given in (11), (12), and (14). They are 

-0.4125 -0.9044 0.1683 -13192 

0.8507 -0.4252 0.3086 -10498 

0.3253 -0.0349 0.9449 430 

-3979 
-15553 
-4251 
















Since a plane can be fixed by any three non-collinear points, we will choose 
vertices 1, 2, and 3. Therefore, from the table in Figure 8, 

(xi. yi, zi) = (0, 0.999,-0.0348), 

(x2, y2* Z2^ “ (0.718, 0.694, -0.034), 

(X3, V3, Z3) = (0.6015, 0.347, 0.719). (16) 

With the aid of (10), (15), and (16) we have the point of intersection given 

by (8): 

^l) = (-0.2004, -0.7835, -0.2141). (17) 

The area enclosed by the window is evaluated using (9): 

= Aj^23 + ^134 = 0.258 + 0.413 

= 0.671 (18) 

As described in article 2.3, we will calculate the areas of the triangles 
whose vertices are (1, 2, I), (2, 3, I), (3, 4, I), and (4, 1, I). They 
are: 

^ Aj^2i ^231 ^341 \ll 

= 0.6725. (19) 

With (18) and (19), we have: 

- a’’ - 0.0042. (20) 

Also using (6), we find that the distance between the intersection 
point (I) and the object (0) is less than that between the pilot (P) 
iid the object (0), i.e.. 


bIO = 

16606.37, 


dPO = 

16607.20. 

(21) 


With the tolerance level e = 0.02 and the table in. figure 5, we conclude 
that the object is visible to the pilot through the window (X) . 



4 . 


CONCLUSION 


We have presented a technique and a computer program to determine whether a pilot 
can see an exterior (point) object if the position and the orientation of the 
^^rcraft, the position of the object, and the geometry of the windows of the 
aircraft are known. The formulation and calculation are extremely simple and 
easy to Implement. The object is treated as a point mass, and may be another 
aircraft, the top of a building, or a mountain peak in the range of sight. 

Pilot look angles to the other references such as runway approach lights 
etc., can also be evaluated. 

Since the speed of the aircraft is finite, and the object most probably 
has a finite dimension, it may be interesting and useful to know thecduratlon 
the pilot can see the object. It might be Interesting as well to determine 
the times and durations the pilot can see the particular area at which he 
expects to land during various traffic patterns. 

The model assumes that the probability of the pilot looking through one 
window is the same as that when looking through another. In practice this is 
not true. The pilot spends more time looking through front windows and less 
through farther ones. This observation may be represented by a Gaussian 
distribution, and included in the calculation of the probability that the 
pilot can see through a particular window. This can be achieved by a 
modification of (Xp, Yp, Zp) at t = t in Equation (8). 


5 . 
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this is The main program 


CVTPST 


This is the main program 


WRITTEN BY J.SaSHIR 
ON machine HONEYWEUL625 


VISION TEST PROGRAM, 

THE MAI'! program calls subroutines INIT AND VISION, 

subroutine init stores the data por the Type or a/c 

CHCSEN iND evaluates ITS FIELD OF VISION, 
subroutine vision determines THE VISIBILITY 
OF A Target fro'* an aircraft whose window 
GEOMETRV IS Known, GIVEN the position coordinates 
of The 'MEhING and Target a/C»S and the yaw, 
pitch A'O bank angles of THt VIEWING A/C, 

input to init is iITypc' gIVEN AS aN ARGUMENT, 

input to SURROUTINE vision is passed through a 
labeled common block • INPUT* .whose Variable 
list is as •- 

vRNWiZ - Runway azimuth is given in radar 
relative FRa'IE of REfFRENCES 
(ZERO DEGREES AT MAG. NORTH), 

VXP,VYP,V?P - POSITION coordinates of the 

viewing A/C<RUNWaY RELATIVE). 

VaZ - VIEWING A/C*s azimuth, GIVEN IN DEGREES 
IN Radar relative coord, system and is 

meaSupED CLCCowISfc From nORTh OR Y.AXlS, 

'VEL - viewing A/C'S elevation angle in degrees; ■■ 
VBA - viewing A/C'S bank ancle in degrefs. 
VXT.VYT.vZT . position coordinates OF THE 
target a/C(RuNWAY RELATIVE). 


OUTPUT IS brought THROUGH A LABELED COMMON BLOCK 

•output* , Whose variable list is as given below, 

ivisibl » 1 means visible, 

» 0 means not visible, 

IVWINOO « window no, aT which target is located. 
_ VPtAMDA « vertical ANGuE OF ThE LOCATION _ 

of the target a/c in the window 



H 
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this is The main procrah 


(UPWARDS IS *VE, downwards -VE) 
VPM'J ■ horizontal ANRlE of the location 
or the target a/c in the wihdow 

(RIGHT IS *VE,UEFT .VE) 


COMMON /INPUT/ VRNWAZ, VXP, VYP,VZP» 

1 VAZi VEL i VBAiVXTi vtt.vzt 

COMMON /OUTPUT/ IViSiBu, IVWpDO.VPLAMOA.ypMU 

read type OF AIRCRAFT DESIRED, 

I’YPE • A CODE NUM3ER WHICH INDICATES 
WHICH A/C IS desired, 

CODES FnR different A/C'S ARE LISTED 
IN SUBROUTINE INIT, 


. . read, ITvPE 


call subroutine INIT TO SET THE A/C TYPE AND 
DETERMI-E ITS FIELD OF VISION, 

call INIT(ITYPC) 

TEST example for FIVE CASES, ‘ " " 

DO 10 I»l,5 

read,vrnwaz,vxp,vyp,vzp, vaz,vel,vba.vxt;vyt,vzt 

call subroutine vision ' 

CALL VISION 

HRITE(6,53)IVISIBL.IVHIND0»VPLAMDA,VPMU 
FCRmat(5x; • IVIS19L « «,I1,AX,'IVWI.ND0 ? 

1 12, AX, »VPLAMDA » • ,F7,2, AX, 'VPMU • *,F7,2,//) 

continue 

STCP 

END 
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16 
17 
19 
19 
2C 
21 
22 
23 
2 * 

25 

26 
27 
23 
29 
33 

31 

32 

33 

34 

35 

36 

37 
39 
39 
*0 






Su^aCUTIVE iNlTdTYPE) 

THIS SuaPOUTiNE STORES THt DATA REQUIRED 
FQR THE TYPE OF AIRCRAFT CHOSEN FOR YISIB- 
-!UITV model, . . .. 

evaluatcn cf the field of vision 

THE field or VISION UTILIZES ALL AVAILABLE 

WINDOWS, aUT EXCLUDES ALL OBSTRUCT 1 CNS , 

THE available field OF VISION IS SUBDIVIDED 
IN-^O NU-BEH Or GEOMETRIC CLOSED FIGURES CaLLEO_ 
windows. 

THfc VER-ICES or THE WINDOWS ARE DETERMINED 
THPOUGH A set cr pats emanating from the PILOT’S 
EYES AN- Each ray IS SPECIFIED IN A SPHERICAL 
COCROlNiTE SvstEM which MEASURES ANGULAR 
DISPLAC=MENTS TPOM FRqm LEVEL REFERENCE AXES 
parallel To The fuselage and parallel to THE 
AIRCRAFT’S WlNGS, 


determination of THE EGUATION OF THE PLANE OF THE WINDOWS 

onlt 3 ••ON collinear points are Required to 

CETERhI''E the fon; of the plane, 

eon AX*ay*CZ*0 DEFINES A PLANE WHERE A,B,C AND D 

ARE constants, THESE CONSTANTS ARfe EVALUATED PY EXPANDING 

the following detfrminant, 

• X Y Z • 

• XI T1 Zl • B 0 

• X2 Y2 Z2 • 

• X3 Y3 Z3 • 


calculation Cf THE AREAS ENCLOSED BY THE WINDOWS. 

eacr window Is divided into NuMBkR OF triangles and 
the area of Each TRUnGLE IS CALCULATED AND THEN 
ADOtO TO EVALUATE THE AREA 0^ THb WINDOW, 

INPUT- 
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label init page . . 2 


43 

44 

45 


17YPI 


- A CODE number which INDICATES 
WHICH A/C IS desired. 


46 


MOONEY. 21 

■ 

1 

SEL 




47 

m 

___ w 

CHEROK£E-1<0 

3 

2 

SEL 




48 

m 

CESSNA. 172 

B 

3 

SEH 




49 

c 

CESSNAa2lO 


4 

SEH 




5C 

m 

w 

PIPER AZTEC 


5 

SEH 




51 

c 

PIPFR.140B 


6 

SEH 



* 

52 


BEECH BaRQN 

U 

7 

TE 




53 

P 

aE»O-COmm(,80E 


9 

TE 




54 

c 

ROEV3i,707 

s 

9 

ME 




55 

• 

DOUGuAS Dee 

- 

10 

me 





56 

57 

58 
. 59 

63 

61 

62 

63 

64 

.. 65 
66 
67 
69 
69 
73 
„71 

72 

73 

74 

75 

76 
_77 

78 

79 
63 
81 

92 

93 

94 


OUTPUT- 

THE OUTPUT VALUES A»E RETURNED IN COMMON BLOCK, 

«;INCE They A«E reouired by subroutines vision 
AND visibl; 

NWIROC - NO. or GEOMETRIC WINDOWS CONSTRUCTED 
OUT or Tm£ available WINDOWS rOR 
Each a/C, 

nrays - Nr or vertices or each window, 

TLaMDa • vertical angle FOR EACH CORNER 
or the window, (IN DEGREES) 
tmu - horizontal angle for Each corner 

or THE window, (IN DEGREES) 


THIS PR-GRAM is written with room for 
expansion to 20 aircrafts, 


COMMON TLAMDA(25,6),TyU<25,6>,NHlND0,XRAy(25,6), 

1 YRAY(25,6),ZRAY(25,6),NRAYS{25), A(4,25), 

2 APEAl(25) 

dimension KW(20) , nr (25,20 > ,Tl< 6, 25,20 ) ,Tm( 6, 25,20 ) , 

1 Vl(l50) , V2(150) , v3(150), V4(i50) , V5(150) , 

2 V6(150) , V7(150) ,V8{150) , V9(150) , VlOdSO), 

3 h1(150),h2(150),h3(150),H4(150),H5(15C), 
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